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Abstract 

Secondary photons from decays of metastable neutralinos can contribute to the 
ultra-high energy cosmic ray flux. The neutralino production rate is too low in ac¬ 
celeration mechanisms to affect the cosmic ray spectrum without emitting enormous 
energy in photons and neutrinos. However, in top-down models with sources not 
concentrated in galactic halos, neutralino decays change the spectrum significantly. 

We estimate the parameters of a model in which photons from neutralino decays 
are responsible for cosmic ray events with energies above 10^*^ eV, and figure out 
distinctive experimental signatures for this model. 

1 . Current experimental data on ultra-high energy (UHE) cosmic ray (CR) spectrum 
are controversial. The results of the AGASA experiment |l[] indicate the absence of the 
GZK feature [^, a cutoff in the spectrum at energies above a few times 10^® eV due to the 
attenuation of high energy protons on cosmic background radiation. However, recently 
published spectrum obtained by the HIRES experiment |P] exhibits this cutoff. Other 
experiments had either low statistics or insufficient precision to clearly support either the 
absence or the presence of the GZK cutoff in the spectrum. 

On the other hand, all experiments (including HIRES) have observed cosmic rays 
with energies as high as (1-^3)- eV . Recently, impressive correlations have been 
found [Q between the arrival directions of UHEGRs registered by the AGASA and Yakutsk 
experiment^ with the positions of gamma-ray loud BL Lac type objects. This implies that 
a signihcant fraction of the extremely energetic particles may originate at cosmological 
distances. Attenuation on the background photons then excludes protons and photons as 
possible candidates for these particles. 

One of the suggested ways to solve the problem is to consider other particles which do 
not attenuate signihcantly on the photonic background. The energy attenuation length 
of photons is of the same order as that of protons. Neutrinos are obvious candidates, but 
neutrino primaries are excluded by the atmospheric shower development [^. However, 

*Based on the talks given by S.T. at Quarks’2002 (the 12th International Seminar on High Energy 
Physics, Novgorod, Russia, June 1-7, 2002) and at SUSY’02 (the 10th International Conference on 
Supersymmetry and Unihcation of Fundamental Interactions, June 17-23, 2002, Hamburg, Germany). 

’^The HIRES stereo data has not been published, nor are we aware of any analysis of the data. 
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UHE neutrinos can scatter off the relic antineutrinos (and vice versa) via the Z-resonance. 
The sites of these so-called Z-bursts serve as secondary sources of photons and nucleons of 
somewhat smaller but still very high energy. If these scatterings take place at a distance 
from the Earth less than the nucleon’s and photon’s energy attenuation length, the Z 
decay products could contribute to the observed UHECR flux [^. The resonant energy 
of a neutrino with mass my is 


F! ~ 

-^res ~ 


4 eV 


m. 


■ 10^^ eV. 


( 1 ) 


Other proposals for non-attenuating particles include rather exotic light supersymmetric 
hadrons and light axion-like particles |^. 

In this talk, we propose an alternative to the Z-burst mechanism which shares a num¬ 
ber of its features but does not require the extraordinarily high particle energies obtained 
from Eq. (®- Like the neutrino, the neutralino can travel for cosmological distances 
unattenuated 0. The lightest superpartners of the Standard Model particles can, in 
some supersymmetric models, decay either to lighter gravitino and non-supersymmetric 
particles (if R parity is conserved), or to non-supersymmetric particles alone with R parity 
violated. These decays can occur more frequently in our cosmological neighbourghood if 
the lifetime of the particle in the laboratory frame is of order but somewhat less than the 
age of the Universe. This gives rise to the super-GZK secondary particles in a way anal¬ 
ogous to the Z burst mechanism. We note that neutralino-induced atmospheric showers 
would be very similar to those induced by neutrinos and hence can be excluded on an 
equal footing to the neutrino events. 

The dominant decay mode of the neutralino is photonic in the models we consider 
here. This means that a signature of this mechanism is the presence of photon-induced 
atmospheric showers. Current experimental data restrict the fraction of photonic showers 
to be less than (28-i-48)% at the energies E < eV [^. However, at higher energies 


the bounds are much weaker, (50 -i- 67)% El. We will see below that the most probable 
implementation of our mechanism is to explain the super-GZK events in the framework 
of a top-down mechanism while relating the events below 10^*^ eV to protons accelerated 
in active galaxies HQ. 

2 . To estimate the required neutralino lifetime and flux, we roughly approximate the 
decay rate of neutralinos as well as the rate of energy loss of photons to the exponentials of 
the distances travelled by particles. We denote the width of decay (neutralino —>• photon 
-|-...) measured in the laboratory frame as T; and the mean energy attenuation length 
of a photon on the cosmic IR and radio background as I ~ 100 Mpc. We suppose that 
the sources are distributed in the Universe with the evolution index m in the comoving 
frame, 

dn{r) = noAnr"^ (1 -|- z{r))'^ dr, r < R, 


where r is the distance from the Earth, z{r) is the corresponding redshift. 

The total UHE photon flux on the Earth, n.y, can be expressed via the total neutralino 
flux from all sources, n^, as 
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777 

-2 

0 

+2 

R, Gpc 

1 

4 

2 

maximal n^/n^ 

0.086 

0.094 

0.068 

T at maximal n^/n^ 

0.5 

6.8 

2.9 


Table 1: R and m are parameters of the source distribution, r is neutralino lifetime in 
the rest frame in units of 10® s • M is the neutralino mass. 


where i?o ~ 4 Gpc is the radius of the Universe, and we calculated the fluxes in the 
laboratory frame. For given m and R, n^/n^ has a broad maximum as a function of F, so 
the hne tuning of F need not be very strong. We present in Table 1 the neutralino lifetimes 
for three sets of values of distribution parameters. Note that the presence of a particle 
with lifetime > 10^ s which decays to photons can affect Big Bang nucleosynthesis due to 
subsequent photodisintegration of light nuclei |^| unless the reheating temperature, and 
hence the particle abundance, are low enough. 

3. Let us turn now to specihc supersymmetric models with metastable neutralino. 
They consist of models with R parity breaking where neutralino LSP can decay to non- 
supersymmetric particles and models with gravitino LSP with conserved R parity (these 
include gauge-mediated supersymmetry breaking (GMSB) |]^ and certain no-scale su¬ 
pergravity models []^). In what follows, we will concentrate on GMSB scenario. 

The lifetime of neutralino-NLSP in the restframe is 


Ibvr^ 


r = 


cos^ 9w ’ 


where 6w is the weak mixing angle, M is the neutralino mass, and F is the scale of 
dynamical supersymmetry breaking. The latter is related to the gravitino mass, 777 . 3 / 2 , as 


F = V3M, 


"^ 3 / 2 ; 


M* = 2.4 ■ 10^® GeV is the reduced Planck mass. We obtain 


F = 2.8- 10 


« GeV 


M 

50 GeV ) 


m,/, = 6,5GeV(^)‘''' 


M 


\ 5/2 
50 GeV / 


( 2 ) 


( 3 ) 


The gravitino is stable due to R parity conservation and its mass is constrained by the 
condition that relic gravitinos do not overdose the Universe |]^. For 7773/2 in the GeV 
range and for reheating temperature low enough to satisfy nucleosynthesis constraints on 
r, the overclosure constraints are satished as well. 

The values (|^), (|^) are on the upper margins for usual gauge mediation but can be 
natural in models of direct gauge mediation. Indeed, let us consider probably the simplest 
complete model of GMSB |M. There, supersymmetry breaking is communicated directly 


from the strongly interacting sector to the MSSM, and 


M 


5 F 

OTT S 


3 


























where s is the messenger scale, and ai(s) is the U{1)y coupling constant taken at the scale 
s. For M 50 GeV and the values of F obtained above using Eq.(|]), this corresponds 
to s ~ 10^'^...10^® GeV depending on the required neutralino lifetime. These values of 
s are within the region allowed for the model of Ref.[jT^ and low enough to suppress 


supergravity contributions to soft masses with respect to GMSB contributions. 

4. The mechanisms responsible for creation of UHE particles can be divided into three 
groups with distinctive observational signatures: 

(1) acceleration in astrophysical sources - arrival directions point back to the sources, 
GZK cutoff is present in the spectrum assuming cosmological distribution of the sources 
and protons or photons as UHE particles; GZK cutoff is absent assuming non-attenuating 
UHE particles. This option seems to be favoured by the data at energies below 10^° eV; 

(2) the decay of metastable relic heavy particles or of short-lived heavy particles orig¬ 
inating in turn from the decay of metastable topological defects which are distributed 
following the Gold dark matter (GDM) density: the sources are concentrated in the ha¬ 
los of galaxies, and the dominant contribution to the observed UHEGR flux comes from 
the halo of the Milky Way (above GZK energy, about 97% for nucleons and photons or 
(15 -i- 30)% for non-attenuating UHE partcles) |]^. Distribution of arrival directions ex¬ 
hibits large-scale anisotropy due to the non-central position of the Sun in the Milky Way 
|r7| . GZK cutoff is absent in the spectrum; 

(3) the decay of short-lived heavy particles originated in turn from the decay of 
metastable topological defects which do not follow the GDM distribution (an example 
of a topological defect which does not follow the GDM density but is distributed more or 
less homogeneously is provided by cosmic necklaces |^): arrival directions of GRs are 
distributed uniformly (unless there are only a few topological defects in the Universe), 
partial GZK cutoff is present for protons or photons; it is absent for non-attenuating 
particles. 

We now consider different mechanisms of neutralino production and check whether the 
mechanisms can produce the required UHEGR flux and not violate other observational 
constraints. We will see that in acceleration mechanisms, option (1), required neutralino 
flux can hardly be produced. 

Indeed, the most probable mechanism of production of neutralinos in astrophysical 
accelerators is in proton-proton collisions. For instance, this could occur in the hot spots 
of active galaxies. All produced supersymmetric particles decay promptly to NLSP. To 
calculate the total neutralino production cross-section, one thus has to sum over all super- 
symmetric species. A collection of expressions for cross-sections can be found in Ref. , 
and approximations for parton distributions can be extrapolated from Ref. pij|. At the 


energies relevant to UHE production, the dominant SUSY production channel is gluon 
fusion. We have checked that the partial cross-section crsusv/o'pp ~ 10“® at these energies, 
where we extrapolated the total pp cross-section from Ref. . If the UHE protons do not 
escape from the source before collision with soft protons (this is the case, for instance, in 
the hot spots of active galaxies), then the total flux of UHE protons in all sources should 
exceed the observed UHEGR flux by a factor of about 10®: 


Hp ^ 


f CrSUSY\ % ir,9 

- —n^ ~ 10 


V ^pp / 
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The observed UHECR energy flux at energies E > 10^° eV is 


E^Jcr{E) 


1 eV cm ^ s ^ 


sr 


-1 


The protonic flux of Jc_r(-E) is excluded for the following reasons. Firstly, the pro¬ 
tons lose their energy by GZK mechanism but do not disappear. Instead, they contribute 
to the CR flux at lower energies [^. The energy flux of protons at sub-GZK energies is 
well measured and is only 10 times larger than the flux at 10^° eV. Secondly, the dominant 
part of the energy flux in pp collisions is released roughly in equal amounts into photons 
and neutrinos - decay products of multiple vr mesons. The photons lose their energy in 
electromagnetic cascades and contribute [Q to the gamma ray background measured by 
EGRET which allows for Jp/JcR ^ 10^ and not 10®. The constraints connected to 
protons and photons can in principle be evaded by very high densities in the sources, so 
that the protons do not leave the source at al^. However, neutrinos cannot be absorbed 
and overshoot the current experimental limits (see Ref. for a recent compilation of 


data) by three orders of magnitude. The only possibility to avoid neutrino production 
is to have enormous densities in the source, ~ 10^® protons/cm“^. Then charged pions, 
which carry about 2/3 of the energy of the products of pp collisions, would interact be¬ 
fore their decay and lose energy in pionic cascades so efficiently that neutrinos would be 
emitted only with low energies. These neutrinos would contribute to a larger atmospheric 
neutrino flux. These densities are hardly possible in realistic astrophysical sources. 

5. We conclude that in the context of the acceleration mechanism, metastable neu- 
tralinos are irrelevant for UHEGRs. On the other hand, in the ’’top-down” mechanisms, 
supersymmetric particles (which promptly decay to neutralino in our case) can carry 
about 40% of the energy of the original heavy particle iE3- Photons from late neu¬ 
tralino decays affect significantly the UHEGR spectrum in the case of homogeneously 
distributed sources (case (3)). The partial GZK cutoff inherent in these mechanisms is 
washed out because neutralino decay probability is higher in our cosmological neighbourg- 
hood. Gurrently, only a limited number of models of the type (3) are marginally consistent 
with EGRET measurements of gamma ray background (see Ref.[^ for examples of 
such models). With metastable neutralinos, EGRET constraints are easily satisfied. The 
mechanism discussed here has the following distinctive signatures in future experiments: 

— a neutralino which does not decay in the detector at future colliders but does not 
constitute the GDM; 

— the absence of positional correlations of GRs with specific astronomical objects at 
energies E > 10^® eV; 

— global isotropy of arrival directions (including absence of galactic anysotropy) at 
E > 102® eV; 

— high fraction of photons at E > 10^® eV. 

In the case (2) of GDM-like distribution of the sources, the dominant part of the UHE¬ 
GRs originate from decays of heavy particles within the Milky Way. Unstable neutralinos 
can affect observable features of GRs in this case only if they decay within the halo, that 
is their lifetime at rest is less than ~ 100 s. 


^ A similar problem appears for the Z burst mechanism if one assumes neutrino origin from pp collisions 
in astrophysical accelerators. In this case, Jp/ JcR ~ 10"'’ is required, and the on-site absorbtion can help. 


though 126 not in the most probable astrophysical sources. 
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